No Stability Switching at Saddle-Node Bifurcations of Solitary Waves in Generalized 

Nonlinear Schrodinger Equations 
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Saddle-node bifurcations arise frequently in solitary waves of diverse physical systems. Previously 
it was believed that solitary waves always undergo stability switching at saddle-node bifurcations, 
just as in finite-dimensional dynamical systems. Here we show that this is not true. For a large 
class of generalized nonlinear Schrodinger equations with real or complex potentials, we prove that 
stability of solitary waves does not switch at saddle-node bifurcations. This analytical result is 
confirmed by numerical examples where both soliton branches are stable at saddle-node bifurcations. 
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Saddle-node bifurcation is a well known phenomenon 
in finite-dimensional dynamical systems In this bi- 
furcation, there are two fixed-point branches on one side 
of the bifurcation point and no fixed points on the other 
side, and the stability of these two fixed-point branches 
switches at the bifurcation point (one branch stable and 
the other branch unstable). In nonlinear partial dif- 
ferential equations (which can be viewed as infinite- 
dimensional dynamical systems), this bifurcation exists 
as well (it is also called fold bifurcation in the litera- 
ture). For instance, solitary waves in nonlinear physi- 
cal systems often exhibit this type of bifurcation. Ex- 
amples include the Boussinesq equations and the fifth- 
order Korteweg-de Vries equation in water waves 
the Swift- Hohenberg equation in pattern formation 5j, 
the nonlinear Schrodinger (NLS) equations with local- 
ized or periodic potentials in nonlinear optics and Bose- 
Einstein condensates [6|-|9(, and many others. Moti- 
vated by stability switching of saddle-node bifurcations 
in finite-dimensional dynamical systems, it is widely be- 
lieved that in nonlinear partial differential equations, sta- 
bility of solitary waves also always switches at saddle- 
node bifurcations (sec 5 8] for examples). This belief 
is very pervasive since no counterexample has been re- 
ported yet. 

In this paper, we show that this belief of universal sta- 
bility switching at saddle-node bifurcations in nonlinear 
partial differential equations is incorrect. Specifically, 
we show that in a wide class of generalized NLS equa- 
tions with real or complex potentials, stability of solitary 
waves actually does not switch at saddle-node bifurca- 
tions. This fact is proved analytically by using general 
conditions of saddle-node bifurcations and eigenvalue- 
bifurcation analysis. It is also verified numerically by 
several examples, where both branches of solitary waves 
are stable at saddle-node bifurcations. 

We consider general nonlinear Schrodinger-type equa- 
tions with arbitrary forms of nonlinearity and external 
potentials in multidimensions, 



iU t + V 2 U + F(\U\ 2 , x)C7 = 0, 



(1) 



where V 2 is the Laplacian in the ./V-dimensional space 
x = (x±,X2, • ■ • ,xn)> an d F{-, •) is a general function 



which contains nonlinearity as well as external poten- 
tials. These equations include the Gross-Pitaevskii equa- 
tion in Bose-Einstein condensates and nonlinear light- 
transmission equations in localized or periodic potentials 
as special cases [Iol - fr2| . Below, we will first focus on the 
case where the function F is real-valued, which applies 
when the system ([T]) is conservative. Extension to the 
non-conservative case of complex functions of F will be 
considered afterwards. 

When the function F is real, Eq. ([T]) admits stationary 
solitary waves of the form 



jut 



it(x), 



(2) 



where u(x) is a localized real function satisfying 

V 2 u- fiu + F(u 2 ,x)u = 0, (3) 

and fj, is a real propagation constant which is a free pa- 
rameter. Under certain conditions, these solitary waves 
undergo saddle-node bifurcations at special values of fi 
. A signature of these bifurcations is that on one side 
of the bifurcation point /iq, there are no solitary wave so- 
lutions; but on the other side of hq, there are two distinct 
solitary-wave branches which merge with each other at 
fi = fiQ. To derive conditions for these bifurcations, we 
introduce the linearization operator of Eq. (J3|) , 



Li 



fx + d u [F(u 2 ,x)u]. 



(4) 



We also introduce the standard inner product of func- 
tions (/, g) = f*(x) g(x)dx, where the superscript 
'*' represents complex conjugation. Our analysis starts 
with the basic observation that, if a bifurcation occurs at 
(i = /iQ, by denoting the corresponding solitary wave and 
the linearization operator as 



u (x) = u(x; ^ ), L w = Li 



M=M0, U=U i 



(5) 



then the linear operator Liq should have a discrete zero 
eigenvalue. This is a necessary condition for all types of 
bifurcations. To derive sufficient conditions for saddle- 
node bifurcations, let us assume that this zero eigenvalue 
of Liq is simple, which is the case for generic bifurcations 
in one spatial dimension as well as for many bifurcations 



2 



in higher spatial dimensions. Under this assumption, we 
denote the unique discrete (localized) eigenfunction of 
I/io at the zero eigenvalue as ip(x.), i.e., 



Lxtfl> = 0. 



(6) 



Since lao is a real operator, we can normalize the eigen- 
function ip to be a real function. We also denote 

G(u;x) = F(u 2 ;x)u, G fc (x) = d*G| u= „ , (7) 

where fc = 1, 2, 3, • • • . Then the sufficient condition for 
saddle-node bifurcations of solitary waves is given by the 
following theorem. 

Theorem 1 Under the above assumption and nota- 
tions, if (ito, "0) 7^ and (G 2 , 4> 3 ) 0, then a saddle-node 
bifurcation of solitary waves occurs at u = uo in Eq. ([3]) . 

Proof. Solitary waves which exist near \x — fiQ admit 
the following perturbation series expansions 



(8) 



Inserting this expansion into Eq. (|3]), we get the following 
equations for Uk at order (n — u ) fe ^ 2 > k = 0, 1, 2, . . . : 



V 2 u - u ""o + F(ul,^.)u 
L w m = 0, 

L w u 2 =u - G 2 u\/2\ , 



(9) 
(10) 
(11) 



and so on. Eq. (j9j for uo is satisfied automatically since 
Mo is a solitary wave at fj, = /io- The ui solution to Eq. 
([TO"! is found from © as 



(12) 



where &i is a constant. The it 2 function satisfies the lin- 
ear inhomogeneous equation (jlip . Due to the Fredholm 
Alternative Theorem and the fact that L10 is self-adjoint, 
Eq. (ITT|) admits a localized solution for u 2 if and only 
if the homogeneous solution if> is orthogonal to the inho- 
mogeneous term, i.e., 



(iP,u -G 2 ul/2) = Q. 



(13) 



Inserting the solution (|12j) into this orthogonality condi- 
tion and recalling the conditions in Theorem 1, we find 
that 



' 2(m ,V>) 



(14) 



we see that these perturbation-series solutions (|15|) give 
two real- valued (legitimate) solitary waves when \x > uo, 
but these solitary waves do not exist when /1 < «o- On 
the other hand, if (uo,tp) and {G 2 ,i/j 3 ) have the opposite 
sign, bf are purely imaginary. In this case, the perturba- 
tion series (|15p give two real- valued solitary waves when 
u < uo but not when u > wo- 

The above perturbation calculations can be contin- 
ued to higher orders. We can show that the two real 
solitary- wave solutions (|15[) . which exist on only one side 
of /1 = uo, can be constructed to all orders of (/1 — uo) 1 ^ 2 - 
In addition, these two solitary waves (x; /i) merge with 
each other when fj, — > (1q. We can also show that ex- 
cept these two solitary- wave branches, there are no other 
solitary- wave solutions near the bifurcation point. Thus 
a saddle-node bifurcation occurs at fi = /xq. This com- 
pletes the proof of Theorem 1. 

Stability properties of solitary waves near saddle-node 
bifurcations is an important issue. In finite-dimensional 
dynamical systems, the stability of fixed points always 
switches at saddle-node bifurcations, and this switch- 
ing is caused by a linear-stability eigenvalue of the fixed 
points crossing zero along the real axis [l|. For solitary 
waves in nonlinear partial differential equations (which 
can be viewed as fixed points in infinite-dimensional dy- 
namical systems), it is widely believed that their stabil- 
ity also always switches at saddle-node bifurcations. We 
find that this belief is incorrect. Below, we show that for 
solitary waves in the generalized NLS equations (flj, 
there are no linear-stability eigenvalues crossing zero at 
a saddle- node bifurcation point, thus stability-switching 
does not occur. 

To study the linear stability of solitary waves © in 
Eq. ([1]), we perturb them as [12| 

U(x,t) = e iflt {u(x) + [v(x)+w(x)]e xt + 

[v*(x) -U!*(x)]e A **} , (16) 

where ii,i»<Cl are normal-mode perturbations, and A is 
the mode's eigenvalue. Inserting this perturbed solution 
into (fT]) and linearizing, we obtain the following linear- 
stability eigenvalue problem 



where 



L 
Li 



$ = 



(17) 



(18) 



Thus, we get two b\ values b 1 which are opposite of each 
other. Inserting the corresponding u\ solutions (|12[) into 
©, we then get two perturbation-series solutions of soli- 
tary waves m(x; u) as 

m ± (x;u) =M (x) + &f(/i-u ) 1/2 V'(x) + 0(m-"o). (15) 

If (u a ,ip) and (G 2 ,ip 3 ) have the same sign, then bf are 
real. Recalling that mq(x) and , 0(x) are both real as well, 



L 



M + ^(m 2 ,x), 



(19) 



and L\ has been given in Eq. (QJ. 

At a saddle-node bifurcation point /i — /iq, we denote 

-^00 = -^0^1=^0, u=u , £-0 = -£^1=^0, ti=tto- (20) 

Then in view of Eq. © , we have 

LqqUq = 0, (21) 
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thus zero is a discrete eigenvalue of Loo- From this equa- 
tion as well as (JSJ), we have 





»o 



= C 



= 0, 



(22) 



thus zero is also a discrete eigenvalue of Cq. 

On the bifurcation of the zero eigenvalue in Co when \i 
moves away from ^o ; we have the following main result. 

Theorem 2 Assuming that zero is a simple discrete 
eigenvalue of Loo and Lio, then at a saddle-node bifur- 
cation point fj>o i n0 eigenvalues of the linear-stability op- 
erator C cross zero, thus no stability switching occurs. 

Proof. The idea of the proof is to show that, when \i 
moves away from ^o > the algebraic multiplicity of the zero 
eigenvalue in C does not decrease, thus the zero eigen- 
value in C cannot bifurcate out to nonzero. 

At the saddle-node bifurcation point fi = /xq, (0,uo) t 
and 0) T are two linearly independent eigenfunctions 
of the zero eigenvalue in Cq in view of Eq. (|22l) . Here 
the superscript 'T represents the transpose of a vector. 
Under the assumption of Theorem 2, zero is a simple 
discrete eigenvalue of Loo and Lio- Thus it is easy to 
see that Cq does not admit any additional eigenfunctions 
at the zero eigenvalue, which means that the geometric 
multiplicity of the zero eigenvalue in Co is two. To de- 
termine the algebraic multiplicity of the zero eigenvalue 
in Co, we need to examine the number of generalized 
eigenfunctions of this zero eigenvalue. The lowest-order 
generalized eigenfunction (/i,<?i) T to the eigenfunction 
(0, uq) t of this zero eigenvalue satisfies the equation 



Co 



so the equation for /i is 



"A" 




' o " 


3i 




u 



Lio/i = u . 



(23) 



(24) 



From Eq. ([6]), we see that this inhomogeneous equa- 
tion has a homogeneous localized solution ip. In addi- 
tion, from conditions of saddle-node bifurcations in The- 
orem 1, {uo,ip) ^ 0. Furthermore, Lio is a self-adjoint 
operator. Thus, from the Fredholm Alternative Theo- 
rem, the inhomogeneous equation (I24[) does not admit 
any localized solution, which means that the eigenfunc- 
tion (0,uo) T of the zero eigenvalue in Co does not have 
any generalized eigenfunctions. Similarly, we can show 
that the eigenfunction (-0, 0) T of the zero eigenvalue in 
Co does not have any generalized eigenfunctions either. 
Hence the algebraic multiplicity of the zero eigenvalue in 
Cq is equal to its geometric multiplicity and is two. 

Away from the bifurcation point (i.e., /i ^ /iq), C al- 
ways has a zero eigenmode 



in view of Eq. (JSJ. In addition, by differentiating Eq. ([3]) 
with respect to [i, we also get 



C 







' " 







u 



(26) 



= 



(25) 



thus (u M ,0) T is a generalized eigenfunction of the zero 
eigenvalue in C. This means that the algebraic multi- 
plicity of the zero eigenvalue in C is at least two when 
M ^ Mo- 

If nonzero eigenvalues bifurcate out from the zero 
eigenvalue in C, the algebraic multiplicity of this zero 
eigenvalue must decrease. Our results above show that, 
when fj, moves away from fj,o, the algebraic multiplicity 
of the zero eigenvalue in C does not decrease, thus there 
cannot be nonzero eigenvalues of C bifurcating out from 
zero. Consequently, no eigenvalues of C cross zero at the 
saddle-node bifurcation point, thus no stability switching 
occurs. This completes the proof of Theorem 2. 

Now we discuss the case when Eq. (fTJ) is nonconserva- 
tive, i.e., the function L(-, •) in {!]) is complex-valued. 
In this case, if F admits parity-time (PT) symmetry 
F*(\U\ 2 ,x) = L(|[/| 2 ,-x), then solitary waves © can 
still exist over a continuous range of real fi values [l3l |. 
and saddle-node bifurcations can also occur (see later 
text). By slightly modifying the analysis above, we can 
show that there is no stability switching at saddle-node 
bifurcations in these nonconservative systems either. 

Next we use two examples to confirm the above ana- 
lytical findings. 

Example 1. Consider Eq. |T]) with a symmetric 
double- well potential and cubic-quintic nonlinearity, i.e., 

iU t + U xx - V(x)U + \U\ 2 U - \U\ A U = 0, (27) 

where the double-well potential 

V[x) = -3 [sech 2 (a; + 1.5) + sech 2 (a; - 1.5)] (28) 

is shown in Fig. (Ha), and the quintic nonlinearity has the 
opposite sign of the cubic nonlinearity. Solitary waves in 
this conservative system are of the form ©, where u(x) 
is real. We have computed these solitary waves by the 
Newton-conjugate-gradient method [Hj], and their power 
curve is plotted in Fig. [TJb). Here the soliton power P 
is defined as \u\ 2 dx. It is seen that a saddle- node 
bifurcation occurs at no ~ 2.16. Two solitary waves on 
the lower and upper branches near this bifurcation point 
are displayed in Fig. [TJc,d). To determine the linear 
stability of these solitary waves, we have computed their 
whole linear-stability spectra by the Fourier collocation 
method [l2j ■ These spectra for the two solitary waves in 
Fig. [Uc,d) are shown in Fig. Q](e,f) respectively. It is 
seen that none of the spectra contains unstable eigenval- 
ues, indicating that these solitary waves on both lower 
and upper branches arc linearly stable. We have also 
performed this spectrum computation for other solitary 
waves on the power curve of Fig. [T^b), and found that 
they are all linearly stable. Thus there is no stability 
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FIG. 1: (a) potential (I28[) in Example 1; (b) power curve of 
solitons; (c,d) soliton profiles at points marked by the same 
letters in (b); (e,f) stability spectra of solitons in (c,d). 

switching at the saddle-node bifurcation point, in agree- 
ment with our analytical result. 

Example 2. We still consider Eq. ([27| but now with 
a complex PT-symmetric potential 

V(x) = -3[sech 2 (a; + 1.5)+sech 2 (a;-1.5)] 

+0.25i [sech 2 (a; + 1.5) - sech 2 (x - 1.5)] , (29) 

see Fig. [Ha). This nonconservative system still admits 
solitary waves ([2]) for continuous real ranges of /i, but 
u{x) is complex- valued now. We have numerically ob- 
tained a family of these solitons, whose power curve is 
plotted in Fig. HJb). Again a saddle- node bifurcation 
can be seen at /iq rs 2.02. For solitary waves on the lower 
and upper branches near this bifurcation point (see Fig. 



[2jc,d)), their stability spectra lie entirely on the imagi- 
nary axis (see Fig. [2][e,f)), indicating that they are all 
linearly stable. Hence no stability switching occurs at 
saddle-node bifurcations in this nonconservative system 
either. 



o 

V 

-2 





(a) 


1 




(C) 


2 




1 










#*\ 












1ft 


r 


uo 


1 

* / 









i 










\l 






















-2 




\ 




-5 


5 


-1 


-5 


5 


- 







• 


X 






X 














(b) 


1 




(d) 




® 


! 


(f) 
























uo 


J? 













\ 




























-2 




i 




1.9 1.95 


2 2.05 


-5 


5 


- 







■ 








X 













FIG. 2: (a) PT potential (gHJ) in Example 2; (b) power curve 
of solitons; (c,d) soliton profiles at points 'c,d' in (b); (e,f) 
stability spectra of solitons in (c,d). In (a,c,d), solid line is 
the real part, and dashed line is the imaginary part. 

In summary, we have shown that for solitary waves 
in the generalized NLS equations ([1} with real or com- 
plex potentials, stability does not switch at saddle-node 
bifurcations. This disproves a wide-spread belief that 
such stability switching should always occur in nonlinear 
partial differential equations. Since the generalized NLS 
equations ([1]) arise frequently in nonlinear optics, Bosc- 
Einstein condensates and other physical disciplines, our 
finding could have broad impact. 
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